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I. INTRODUCTION 

The quarkonia, especially the bottomonium bb states, are approximately non-relativistic 
systems since they do not contain intrinsically relativistic light quarks. Hence, these states 
are the best candidates to examine the hadronic dynamics and investigate both perturbative 
and non-perturbative characteristic of QCD. In the past, mainly theoretical calculations 
on the properties of these states had been made using potential model or its extensions 
like the Coulomb gauge model (see for instance [iFs} and references therein). In jjj], both 
potential model and QCD sum rule approach have been applied to extract the ground- 
state decay constant of mesons containing heavy b quark. It is stated that the QCD sum 
rule technique gives more reliable and accurate determination of bound-state characteristics 
compared to the potential models by tunning the continuum threshold parameter. The QCD 
sum rule approach [8| is one of the most powerful and applicable tools to hadron physics. 
This model has been widely applied to investigate the spectroscopy of hadrons and their 
electromagnetic, weak and strong decays. The obtained results have very good consistencies 
with the experimental data to date within the typical (10-20)% error bars of the technique. 

The present work is dedicated to investigation of the semileptonic transition of scalar 
P wave bottomonium xm(lP) meson with quantum numbers I G (J PC ) = + (0 ++ ) into the 
pseudoscalar B c meson. The x&o(l-P) state has been observed first in radiative decay of the 
T(2S) 0] and recently has been confirmed by ATLAS Collaboration [lo| together with the 
higher x b (2P) and %b(3P) states. In the latter, these quarkonia states have been produced 
in proton-proton collisions at the Large Hadron Collider (LHC) at y/s = 7 TeV and through 



their radiative decays to T(1S,2S) with T — > fi + [i~. Our previous theoretical results [U| 
on the mass of these states done both in vacuum and finite temperature QCD are in good 
agreement with the experimental results [loj]. Note that, we also have applied the QCD sum 
rules approach both in vacuum and finite temperature to investigate the spectroscopy of the 
pseudoscalar, vector and tensor quarkonia in [l2-14]. As we know the masses and decay 



constants of the quarkonia, it is possible to investigate their electromagnetic, weak (leptonic- 
semileptonic) and strong decays. Considering such decay channels can help us obtain more 
information about the nature of the scalar xwO-P) meson as well as perturbative and non- 
perturbative aspects of QCD. 

The layout of this article is as follows. In the next section, we derive the QCD sum 
rules for the form factors appearing in the amplitude of the semileptonic decay channel 
under consideration. To do so, we take into account the two-gluon condensates as the 
non-perturbative contributions to the correlation function. Section III is devoted to our 
numerical analysis of the obtained form factors and their behavior in terms of the transferred 
momentum squared. In this section, we also numerically estimate the decay width of the 
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semileptonic xw B c tv decay mode. The last section encompasses our concluding remarks. 

II. QCD SUM RULES FOR TRANSITION FORM FACTORS OF X bo -> B c £u 

The hadronic event under consideration can be described in terms of quark degrees of 
freedom by the process b — > civ at tree-level, whose effective Hamiltonian can be written as: 

G 

Heff = -^ V cb v 7m(! - c 7 M (1 - 75)6, (1) 

where Gp is the Fermi weak coupling constant and is an element of the Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix. The transition amplitude is obtained via 

M = (B c (p')\H eff \xw(p)), (2) 

or 

M = ^V cb v^(l - l5 )l{B c (p') I c 7m (1 - 75)^ I Xw(p)). (3) 

To proceed, we need to know the transition matrix element (B c (p') | c 7m (1 — 75)6 | Xbo(p)) 
whose vector part do not contribute due to parity considerations, i.e., 

{B c (j/)\cry lt b\ X n(p))=0. (4) 

The axial-vector part of transition matrix element can be parameterized in terms of form 
factors as 

(B c (p') I cj^b I X bo(p)) = fi(q 2 )P, + f2(q 2 )q„ (5) 

where fi(q 2 ) and f2(q 2 ) are transition form factors; and P M = (p + p')^ and q^ = (p — p%. 

Our main goal in the present section is to calculate the transition form factors applying 
the QCD sum rules technique. The starting point is to consider the following tree-point 
correlation function as the main ingredient of the model: 

U, = i 2 J d 4 x J d 4 ye-^'y(0 I T{ JbM < ;F (0) 4»} I 0), (6) 

where T is the time ordering product, Jb c (u) = cj 5 b and J Xb0 (x) = bUb are the interpolating 
currents of the B c and xm mesons, respectively; and J^(0) = cj^b and (0) = ^7^756 are 
the vector and axial- vector parts of the transition current. Following the general idea in the 
QCD sum rules technique, we calculate this correlation function once in terms of hadronic 
degrees of freedom called physical or phenomenological side and the second in terms of 
QCD degrees of freedom (quarks and gluons and their interaction with QCD vacuum) called 
the QCD side. The latter is done in the deep Euclidean region by the help of operator 
product expansion (OPE). These two representations are then matched together, using the 
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quark-hadron duality assumption, through a double dispersion relation to obtain the QCD 
sum rules for the form factors. As we deal with the ground states in this approach, we 
shall separate the ground state from the higher states and continuum. This is done by 
two mathematical operations called Borel transformation and continuum subtraction. Such 
transformations bring some auxiliary parameters namely two Borel mass parameters and 
two continuum thresholds for which we will find their working regions in the next section. 

The phenomenological side of the correlation function is obtained inserting two complete 
sets of intermediate states with the same quantum numbers as the interpolating currents 
J Bc and J XbQ . As a result, we obtain 

^phys = (0 I JbM I I #(0) I Xm(p))(Xm(p) I 4J°) I °) m 

" ~ (p' 2 - ™\ C W - KJ ' 

where • • • represents the contributions coming from higher states and continuum. Besides the 
transition matrix elements defined previously, the matrix elements of interpolating current 
between the vacuum and hadronic states are parameterized in terms of the leptonic decay 
constants, i.e., 

(0 I J Bc (0) I B c {p')) = i-^gs (xw(p) I 4J°) I °) = -im Xw f Xb0 . (8) 

Putting all expressions together, the final version of the phenomenological side of the corre- 
lation function is obtained as 



fxm m xm fB c m 2 Bc 



Y[ PHYS (p 2 p' 2 \ — t XbO ' ' "XW jr>c"»B c 

M ' (p /2 — m B J(p 2 — iri^ ) fnb + m. 



fi(q 2 )P, + f2(q 2 )q, 



+ ..., 

(9) 



where we will choose the structures P M and g M , to evaluate the form factors f\ and / 2 , 
respectively. 

At QCD side, the correlation function is calculated in deep Euclidean region by the help 
of the OPE. For this aim, we write the coefficient of each structure in correlation function 
as a sum of a perturbative (diagram a in figure 1) and a non-perturbative (diagrams b, c, 
d, e, f and g in figure 1) parts as follows: 



U QCD = (jjpert + n „o„ pert)p ^ + ^pert + ^nonpert^ ^ 

where, the nf er< functions are written in terms of double dispersion integrals in the following 

way: 



Ilf 6r * = - ( ds I ds'-. — — S ' S . '- q — -r- + subtraction terms, (11) 

(2%) 2 J J (s — p 2 )(s' — p' 2 ) 
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FIG. 1. Feynman diagrams contributing to the correlation function for the Xbo ~^ B c tU decay: (a) 
the bare loop and (b, c, d, e, f, g) two-gluon condensate diagrams. 

where, pi(s,s',q 2 ) are the spectral densities with i = 1 or 2. Applying the usual Feyn- 
man integral technique to the bare loop diagram, the spectral densities are calculated 
via Cutkosky rules, i.e., by replacing the quark propagators with Dirac delta functions: 



-2ir5(p 2 



m 



implying that all quarks are real. After some calculations, the 



spectral densities are obtained as follows: 

2N c I (s, s', q 2 ) mb(m c — 3m;,) — A(h + s) — B{h + s' 
2NJ (s, s', q 2 ) \m b (m b + m c ) - A(h - s) + B(h - s') 



Pi{s,s ,q ) 



p2{s,s,q 



(12) 



where 



I (s,s',q 2 



4AV2, 



s,s',q 2 ) 



A 

B 
h 



A(s,s',g 2 ) 
1 

\(s,s',q 2 ) 
2m b {m b — m 



( 2 2 \ If i / 

(m b — m c )u + us 



2(mi 



m c )s + su 



(13) 



here also X(s, s', q 2 ) = s 2 + s' 2 + q 4 — 2ss' — 2sq 2 — 2s' q 2 , u = q 2 + s — s' , u' = q 2 — s + s' , 
u" = q 2 — s — s' and A^ c = 3 is the number of colors. The integration region for the 
perturbative contribution in Eq. ffTT]) (bare loop diagram) is determined requiring that the 
arguments of the three 5 functions vanish, simultaneously. Therefore, the physical region in 
the s and s' plane is described by the following inequality: 

2s[m 2 b - m 2 c + s' + *f ] 



-l<f{s,s') 



\ l / 2 (ml s,m 2 )X l / 2 {s, s\q 2 ) 



< +1. 



(14) 
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For the non-perturbative part, we take into account the two-gluon condensate diagrams (b, 
c, d, e, f, g) in figure 1. Here we should mention that we deal with the heavy quarks in the 
present work and the heavy quarks' condensates are suppressed by inverse of their masses, so 
we can ignore them safely. After lengthy calculations for the two-gluon condensates diagrams 
(b, c, d) correspond to the diagrams with two gluon lines coming out from different quark 
lines, we get 



r l rl-x 1 f (2)1 (2)2 q4 



where 



and 



V = m c x — m b r + p xr + pyf + xy(p + p — q ), 



(15) 



(16) 



Q\ = ^^9mlm c y 2 v 2 (-l + 2y) - 18mly 2 v 3 - 3m\xv 3m 2 c (r 2 x + 2rxy + 2wy 2 ) 
+ q 2 y [ 3r 3 x + 3r 2 ty + ry 2 (-U + 15a;) + y 3 (-22 + 18x + 9y) 



+ 3m\m c xv 



+ ni\ ) m c q 2 x 2 yv 



x(Q- I9x + 13x 2 - 15y + 22xy + 12y 2 ) + q 2 y(^13r 2 x 2 - 6r 2 x 
+ 3ry - Ihryx + 22ryx 2 + y 2 (l5 - 22y - 26x + 12x 2 + 12xy + 12y 2 ) 

\x (22 - 49x + 39x 2 - 62y + QQxy + 3Qy 2 ) + q 2 y |^39x 4 - 88x 3 
+ QQx 3 y + 3fy{3- 7y + Qy 2 ) +x 2 (jl - 128y + 3Qy 2 ) + xz(ll - 33y + 18y 2 )^j 
+ 3q 4 x 3 y 2 m 2 c |/ 2 + Qf 3 y + x 3 (-ll + 3x + 9y) + /x(lO - 17a: + 18xy - 26y + 18y 2 )j 
+ q 2 y(^x 3 (31 + 3x 2 - Ux + 9xy - 50y + 21y 2 ) + /x( - 20 + 5% - 51y 2 + 15y 3 



+ 3m\q 2 x 2 y 



q 2 y\Ar 3 - 2r 3 x + l6r 2 y 



+ 36x - 59xy + 27xy 2 } + / 2 (-4 + 3yz 2 ) 
- 8r 2 xy + 3ry 2 (8 -7x + 2x 2 ) + y 3 (l7 -by- 27 x + I2x 2 + 6xy)^ + m 2 ^6xy - 4x 
+ hx 2 - 2x 3 - hx 2 y + 6a; V + f(8y- 8xy + 12xy 2 + Qy 2 z - l^j j, 



0? = — < 9m\m c v 
1 96| 6 



rx 2 (-7 + 13a;) + 2a; 2 (-8 + llx)y - 2y 2 (3/ + 5x - 6x 2 ) + 12y 3 (/ + x) 



9m b v 



3r 2 x 2 + Qrx 2 y + y 2 (j - 14a; + 15x 2 - 16y + I8xy + 9y 2 ) 



+ 3m fe x 
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x 



ey I A r\ Q\ / 

3m\Qy — xr + 6y r — 2xry + \2ry J + q y\8r x — llr x — Ylry + hArxy 
38ryx 2 - y 2 (A3 + 153a: 2 - 48a: 3 + 5Aty + 23y 2 - 158a; - A8xyt - A8xy 2 ^j 

9m 2 a;(4 + I3rx- I3y + 22xy + 12y 2 ) + g 2 y( 312x 4 + x 3 (-535 + 52%) 



+ mbm c xv 



+ a: 2 (313 - 683y + 288?/ 2 ) + 3y( - % + 53y - 88y 2 + 48y 3 ) + 6a; ( - 9 + AOy 
- 58y 2 + 24t/ 3 ) 



+ 3q z x A y 



: |V 2 (3-29a; + 24a; 2 ) - y (21 - x(l51 - 206a: + 72x 2 )^J 

+ 78y + 2xy(-m + 72x) + 36y 2 (-3 + Ax) + A8y 3 ^j + q 2 y (^A5x 5 + a; 4 (-182 + 135y) 
+ a: 3 (306 - 584y + 315y 2 ) + x 2 ( - 253 + 8A9y - lOOAy 2 + 405y 3 ) + 3/(4 - 30y + 69y 2 
- 59y 3 + I5y 4 ) + x(96 - 500y + 965y 2 - 792y 3 + 225y 4 
1 



1 96 \ 6 

+ 9m 2 x 

+ Aq 2 yx 



rV - 2rx 2 y + r(l- 7xy 2 + 8xV) + 2y 3 (l + 8rx) + y 4 (-l + 8a:) 
r 2 a:(-3 + Ax) + 2rxy(-5 + 6x) + 2y 2 (2 - 13x + 12a: 2 ) + 12ty 3 + 8y 
r 2 x(l7 - 64a: + 45a; 2 ) + y(21 + x(-154 + 413a: - 418a: 2 + 135a: 3 )) - 120y 



+ y 2 x(A97 - 718x + 315a; 2 ) + 3y 3 (76 - 198a; + 135a; 2 ) + 3?/ 4 (-58 + 75a;) + A5y z 



+ 3mbm c v 



52a: 4 + 2Axfy 2 + 12/^ 2 (— 1 + 2y) + a: 3 (-61 + 88y) + 3a: 2 (7 - 19y + 16y 2 ) 
6 4 = j^^>x 5 + y 2 f 2 (-A + by) + a; 4 (-14 + 15y) + fxy 2 (-19 + 25y) + x 3 (l3 - 28y + 3hy 2 ) 



x; 



; (4-13y + 48y 2 -45y 3 )|. 



;i7) 



Here v = — 1 + x + y, r = — 1 + x, t — — 1 + 2a:, w — 1 + x, f = — 1 + y and z = —2 + y. 
The explicit expressions for B^' 2,3 ' 4 correspond to the structure g M are too long, hence we do 
not present them here. In a similar way, we calculate the contributions of the diagrams (e, 
f, g) correspond to two gluon lines coming out from the same quark line. Because of their 
very lengthy expressions, we do not also depict their explicit form here, but we will take 
into account their contributions in our numerical results. 

The next step is to equate the coefficients of selected structures from both sides in order 
to get sum rules for the form factors. After applying double Borel transformations with 
respect to the variables p 2 and p' 2 (p 2 — » M 2 ,p' 2 — > M' 2 ) to suppress the contributions of 
the higher states and continuum, the QCD sum rules for the form factors are obtained as: 
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h 2 (g 2 ) = ( ? + r c) , 1 e<J M2 e- 2 B J M ' 2 { - JL p ds r'° ds > 

fB c m 2 Bc f Xbo m Xbo I (27r) 2 7 4 mf ,/(™ 6 +™c) 2 

x p 1)2 ( S , s', q 2 )6[l - f(s, s')]e- s / M2 e' s '/ M ' 2 + Bn n ™ pert }, (18) 

where the operator B denotes double Borel transformation. Note that to subtract contribu- 
tions of the higher states and continuum, we also apply the quark-hadron duality assumption, 
i.e., 

p hi9her states (s, s') = p OPE (s, s')9(s - s Q )0(s' - s' ). (19) 
We also perform the double Borel transformation as follows: 
• for the perturbative part, we use 

B t - 1 21 r - 1 21 11 e -™?/M 2 e - m |/M- 



[p 2 - m\] m \p' 2 ~ mlY v ' r[m]r[n] (M 2 ) m - 1 (M /2 ) n ~ 1 ' 

(20) 

For the non-perturbative part, first we make the transformation 

I e -/(p' 2 )/M 2 

\P 2 -f(p' 2 )} n T[n](M) n - 1 ' K ' 

to write the terms containing p' 2 in exponential form. Then we apply the following 
rule to transform the (p 2 — > M' 2 ): 

^- aP " 2 = ^ " «)> (22) 

where a is a function of quarks' masses as well as the parameters used in Feynman 
parametrization. 



III. NUMERICAL RESULTS 



In this section, we numerically analyze the related form factors, obtain their fit function in 
terms of q 2 and estimate the decay width of the channel under consideration. In calculations, 
we use the input parameters as presented in table I. 

The sum rules for the form factors denote that they also depend on four auxiliary pa- 
rameters, namely continuum thresholds So and s' and Borel mass parameters M 2 and M' 2 . 
The continuum thresholds are not completely arbitrary but they are in correlation with the 
energy of the excited state in initial and final channels. Considering this point and the 
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Parameters 


Values 


m c 


(1.275 ± 0.015) GeV 


m b 


(4.7 ±0.1) GeV 


m e 


0.00051 GeV 


mn 


0.1056 GeV 


m T 


1.776 GeV 


m XbO 


(9859.44 ± 0.42 ± 0.31) MeV 


m Bc 


(6.277 ± 0.006) GeV 


fB c 


(400 ± 40) MeV 


fxbO 


(175 ± 55) MeV 


G F 


1.17 x 10" 5 GeV~ 2 


Vcb 


(41.2 ± 1.1) x 10~ 3 


(0|i«sG 2 |0) 


(0.012 ± 0.004) GeV 4 



TABLE I. Input parameters used in our calculations 




fact that the result of the physical quantities (form factors) should weakly depend on these 
parameters, we choose the intervals so = (97.7 — 99.2) GeV 2 and s' Q = (40 — 41) GeV 2 
slightly higher than the mass of pole squared of the initial and final mesonic channels for the 
continuum thresholds. The Borel parameters M 2 and M' 2 also are not physical quantities, 
hence the form factors should be independent of them. The reliable regions for the Borel 
parameters M 2 and M' 2 can be determined by requiring that not only the contributions of 
the higher states and continuum are effectively suppressed, but contribution of the operators 
with the higher dimensions are small, i.e., the sum rules for form factors converge. As a 
result of these requirements, the working regions for these parameters are determined to be 
15 GeV 2 < M 2 < 30 GeV 2 and 10 GeV 2 < M' 2 < 20 GeV 2 . The dependence of form factors 
fi and /2 on Borel masses at q 2 = 1 GeV 2 are plotted in figure 2. From this figure, we see 
good stability of the form factors with respect to the variations of the Borel mass parameters 
at their working regions. To see the convergence of the OPE, we compare both perturbative 
and non-perturbative contributions to the form factors in figure 3 at q 2 = 1 GeV 2 and the 
presented Borel windows. From this figure and our numerical calculations, it is found that 
the ratio of non-perturbative contribution to that of perturbative is 0.08 and 0.05 for fi 
and f 2 , respectively. Hence, the non-perturbative contribution constitutes only 7.5% and 
4.8% of the total results respectively for the form factors fi and f 2 - This means that the 
series of sum rules for the form factors are convergent. In the presented Borel windows, 
the contributions of the excited and continuum states are exponentially suppressed. This 
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FIG. 2. Dependence of the form factors fa and fa on Borel mass parameters M 2 and M' 2 at 
g 2 = ! Ge y2 




FIG. 3. Comparison between perturbative and non-perturbative contributions to the form factors 
at q 2 = 1 GeV 2 and chosen Borel windows. The upper (lower) plane belongs to the perturbative 
contribution in fa (fa). 

guarantees the reliability of the sum rules and isolation of the ground state from the excited 
states and continuum. 

Our calculations show that the form factors are truncated at q 2 ~ 9 GeV 2 (see figure 4). 
After this point up to the hi ghe r limit of the q 2 , the sum rules predictions are not reliable 
(for details see for instance |l7i 0]). However, we need their fit functions in the whole 
physical region, m 2 < q 2 < (jn Xw — tub,,) 2 to estimate the decay width of the Xbo -> B c £v 
transition. To extend our results to the full physical region, we search for parameterization 
of the form factors in such a way that in the region < q 2 < 9 GeV 2 , predictions of this 
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parameterization coincide with the sum rules results. The following parametrization adjust 
well the q 2 dependence of the form factors: 

f^q 2 ) = a —^- + — 2 (23) 

(I q , ) (I q , ) 2 

v mi ~ I v mi.. I 

fit fit 

where, the values of the parameters a, b and mfu obtained using M 2 = 25 GeV 2 and 
M' 2 = 15 GeV 2 for the Xbo — ► B c iV channel are given in table II. 





a 


b 


mj tt (GeV 2 ) 


fxiXM -»> B c £v) 


-0.055 


0.062 


21.86 


/ 2 (X60 -> ^c^) 


0.225 


-0.254 


19.79 



TABLE II. Parameters appearing in the fit function of the form factors. 

We depict the dependence of form factors f\ and /2 on q 2 obtained directly from the sum 
rules as well as the fit parametrization at whole physical region in figure HI In the case of 
sum rules predictions, we present the perturbative, non-perturbative and total contributions 
in this figure. 



0.15 k, 

Total 
Perturbative 
Non -Perturbative 

0.10 
0.05 

0.00 

2 4 6 8 10 12 2 4 6 8 10 12 

? 2 (GeV 2 ) ? 2 (GeV 2 ) 




FIG. 4. Dependence of the form factors fi and / 2 on q 2 at M 2 = 25 GeV 2 and M >2 = 15 GeV 2 . 

Having obtained the behavior of the form factors in terms of q 2 at whole physical region, 
we would like to calculate the decay width of the process under consideration. Using the 
amplitude previously discussed, the differential decay width for xto — > B c £u is obtained in 
terms of form factors as: 



dT G 2 F |2xl/2/ 2 2 2> 



O 
to 








1 T) \ 


[2 . 



|/ 1 (g 2 )| 2 (2m| c + 2m|-g 2 ) 



12 



+ 2« - m\)Ite{h(q 2 )r 2 {q 2 )] + \f 2 (q 2 )\ 2 q 2 
+ 2{ml - mDq'Reihiq^ttiq 2 )} + \f 2 (q 2 )\ 2 q A 



f + mf\ 



2q 2 



\fi(q 2 m< c -™ 2 s 



2 \2 



Performing integration over q 2 in Eq. 



in the interval mf < q 2 < (w\ b0 — rriB c 



(24) 



we 



obtain the expression for the total decay width. The numerical values of the decay width at 
different lepton channels are presented in Table III. The errors in the values of the decay rates 





T(GeV) 


Xbo -> B c eV e 


1.46 x 10" 14 


XbO -> BcfJA/ft 


1.45 x 1(T 14 


Xbo -> B c tV t 


0.91 x 10" 14 



TABLE III. Numerical results for decay rate at different lepton channels. 

in table III are due to uncertainties in determination of the working regions for continuum 
thresholds and Borel mass parameters as well as errors of the other input parameters. 



IV. CONCLUSION 



In the present work, we studied the semileptonic xto — > B c £v, (£ = (e, /i, r)) decay channel 
within the framework of the three-point QCD sum rules. In particular, taking into account 
the two-gluon diagrams as non-perturbative contributions, we obtained the QCD sum rules 
for the form factors entered the transition matrix elements. After obtaining the working 
regions for the auxiliary parameters, we found the behavior of the form factors in terms of 
q 2 in whole physical region. The fit function of the form factors were then used to estimate 
the decay rates at different lepton channels. Any measurement on the form factors as well as 
decay rate of the channel under consideration and comparison of the obtained results with 
theoretical predictions in the present study can give valuable information about the internal 
structures of the participating mesons specially nature of the scalar Xbo(lP) state. 
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